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Synthetic magnetism in cold atomic gases opened the doors to many exciting novel physical 
systems and phenomena. Ubiquitous are the methods used for the creation of synthetic magnetic 
fields. They include rapidly rotating Bose-Einstein condensates employing the analogy between 
the Coriolis and the Lorentz force, and laser-atom interactions employing the analogy between the 
Berry phase and the Aharonov-Bohm phase. Interestingly, radiation pressure - being one of the most 
common forces induced by light - has not yet been used for synthetic magnetism. We experimentally 
demonstrate a synthetic Lorentz force, based on the radiation pressure and the Doppler effect, by 
observing the centre-of-mass motion of a cold atomic cloud. The force is perpendicular to the velocity 
of the cold atomic cloud, and zero for the cloud at rest. Our novel concept is straightforward to 
implement in a large volume, for a broad range of velocities, and can be extended to different 
geometries. 


INTRODUCTION 

Experiments on synthetic magnetic/gauge fields for neutral atoms [T1 !TTi| have enabled realizations of the Hall 
effect 0, famous Hamiltonians such as the Harper 0, I| and the Haldane Hamiltonian 0, intriguing topological 
effects [9, 0 , and the observation of synthetic Dirac monopoles 0- There are a few recent reviews on this promising 
theoretical and experimental progress in synthetic magnetic/gauge fields (l^ - ll5| . The first implementation of synthetic 
magnetism was in rapidly rotating Bose-Einstein condensates (BECs), employing the analogy between the Lorentz 
force and the Coriolis for ce 0 , Q . The methods based on laser-atom interaction employ the analogy between the Berry 
phase in atomic systems [161 ] , and the Aharonov-Bohm phase for charged particles |l3| . The first of them was realized 
in the NIST group with spatially dependent Raman optical coupling between internal hyperfine atomic states in bulk 
BECs 0. Methods generating synthetic magnetic fields in optical lattices engineer the complex tunnelling matrix 
elements between lattice sites [5H9(. They include shaking of the optical lattice 0, laser assisted tunnelling in optical 
superlattices realizing staggered synthetic magnetic fields Q, in tilted lattices realizing homogeneous fields 0,0, and 
an all-optical scheme which enables flux rectification in optical superlattices Q ■ Interestingly, radiation pressure has 
not yet been used among the methods for synthetic magnetism, that is, to create the analogue of the Lorentz force. 
Here we experimentally demonstrate the synthetic Lorentz force based on radiation pressure in cold atomic gases. We 
measure the dependence of the transverse radiation pressure force (analogous to the transverse Hall deflection) on the 
velocity of a cold atomic cloud by observing the centre-of-mass (CM) motion. The observed force is perpendicular to 
the velocity, and zero for the atomic cloud at rest. This concept based on radiation pressure, theoretically proposed 
in Ref. is straightforward to implement in a large volume (e.g., volumes 1 mm 3 -l cm 3 are easily accessible) fl8[ . 
is applicable for a broad range of velocities, and can be extended to different geometries. The main reason for the 
absence of radiation pressure from the previously used methods of synthetic magnetism is the associated heating due 
to spontaneous emission. However, this is not an obstacle for atomic gases cooled and trapped in a Magneto-Optical 
Trap (MOT), where our experiments are performed. 


RESULTS 

The basic idea. The idea behind our experiment is to drive two-step two-photon transitions between three atomic 
levels, |0) |1) —► 12), with mutually perpendicular laser beams as illustrated in Fig. [Q Two counterpropagating 

laser beams aligned with the a;-axis drive the |0) —> |1) transition, whereas the |1) —> |2) transition is driven by 
counterpropagating beams aligned on the y-axis. Due to the Doppler effect and the perpendicular configuration of 
the laser beams, both components of the radiation pressure force depend on both components of the atomic velocity: 
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FIG. 1. The scheme used to obtain the synthetic Lorentz force via radiation pressure. It is based on two-step two-photon 
transitions involving three atomic levels, (a) Two counter-propagating laser beams on the axaxis (red arrows) drive the |0) —► |1) 
transition, while the |1) —¥ |2) transition is driven by two counterpropagating beams on the y-axis (blue arrows). The transverse 
radiation pressure force arising from the |1) —> |2) transition F y , depends on the velocity v x as indicated, (b) Four possible 
excitation pathways and a sketch of the detuning values: P(—¥,\.) denotes absorption of a |0) —> |1) photon going towards the 
positive ^-direction, followed by absorption of a |1) —> \2) photon travelling towards the negative {/-direction, and so on. The 
total detuning value for P(—>,4-) and P(<— , j") i s A, whereas it is much larger in magnitude for P{— »,t) an d P(<— ,4-); the two 
latter pathways are thus negligible in this configuration, (c) The Doppler shift provides F y as a function of v x as sketched in 
(a): pathway P(— 1,4_) for an atom with positive velocity v x = +A/k x ( v y = 0) is on resonance, providing negative F y . Likewise, 
P(<— j't') is on resonance for an atom with negative velocity v x = —A/k x , providing positive F y . See text for details. 


F x = F x (v x ,Vy) and F y = F y (v x ,v y ). This gives us the opportunity to design the detuning values of our lasers such 
that F y is positive/negative for atoms with negative/positive velocity component v x , and that the total force is zero 
for an atom at rest: F(v = 0) = 0. These are the characteristics of the synthetic Lorentz force that we experimentally 
demonstrate. 

The design of the detuning values of the lasers is crucial in obtaining the desired result. The beams driving the 
first step of the transition |0) |1) are detuned by the same magnitude, but with the opposite sign. The one towards 

the positive ^-direction is red-detuned by = —5\ < 0, while the other is blue-detuned by 5<- = 5 1 > 0. Their 
intensities are equal. Thus, if just these two lasers were present, the net force on atoms (of any velocity) would be 
zero. However, the population of level |1) would depend on the velocity v x , which implies that the rate of transitions 
|1) —> |2) giving the transverse force will depend on v x . The detuning values of the beams driving the second step of 
the transition, |1) — > \2) are denoted by <5q. and <5^, for the beam propagating in the positive and negative y direction, 
respectively. For now, let us set these values such that i5_> + = 5<- + c)-|- = A > 0, as indicated in Fig. [0(b). 

The two-step two-photon transitions, where absorption of a |1) —> |2) photon follows absorption of a |0) —► |1) 
photon with perpendicular momentum, yield the synthetic Lorentz force via momentum transfer from photons to 
atoms. Given the fact that we have two counterpropagating beams for each transition, we have four excitation 
pathways for the two-step two-photon transition, denoted by P{— I,f), P(t—,t)> P(~' K-I)j and P(<—,1), see Fig. [Q 
The arrows correspond to the direction of the photon’s momentum, for example, P(—>,J,) denotes the pathway where 
absorption of a photon travelling in the +x direction is followed by absorption of a photon in the — y direction and so 
on. Since the detuning magnitudes of the first step are identical for all pathways (|<5_>| = &<_), the relevant quantity 
is the total detuning for the two-step two-photon transitions. It is important to note that the total detuning A for 
pathways P(— 1,4-) and P(-s—, t), is much smaller in magnitude than the detuning values of P(— j*) and P(y- ,-J.). The 
last two are thus negligible in the setup of Fig. [TJ 

To understand the origin of the synthetic Lorentz force, we take into account the Doppler shift. For an atom moving 
along the cc-axis with velocity v x ( v y = 0), the P(<—, j') pathway is detuned by A + k x v x , i.e., it is on resonance when 
v x = —A /k x < 0. Because photons from the second step of P(t—,t) impart momentum towards the positive y- 
direction, there will be a positive force F y for atoms with negative v x . In the same fashion, the P(—>,4-) pathway will 
be on resonance when v x = +A/k x > 0 yielding negative F y for positive v x . 

The experiment. In the experiment, we use cold 87 Rb atoms. For the |0) —► |1) transition we use the D2 
transition in 87 Rb: | 5Si/ 2 \F = 2) —>• |5Pw?; F' = 3) at 780 nm [l§|. The easily accessible transition, |5P3/2;P / = 
3) —> \5D 5 /2',F" = 4) at 776 nm [20|, is used for the second step |1) —> |2). The linewidths of the states are 
27 t x 6.1 MHz for the |5P3/2;P' = 3) state [l9|, and 27r x 0.66 MHz for the \5D 5 /2',F" = 4) state [20| . 
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The 87 Rb atoms are cooled and trapped in a standard glass vapour cell magneto-optical trap (MOT), arranged in a 
retro-reflected configuration. In this configuration, three orthogonal retro-reflected beams are used to create the total 
of six beams needed for the MOT. A pair of anti-Helmholtz coils provides a quadrupole magnetic field, which together 
with the laser beams creates a trapping potential (for example, see Ref. [18| ). Fluorescence imaging of the cloud is 
performed with a camera aligned along the z-axis. In typical experimental conditions we obtain a cloud of 0.4 mm in 
diameter, which contains about 10 s atoms of 87 Rb, at a temperature of 50 yiK (for details of the experimental setup 
see Methods). The four beams implementing the synthetic Lorentz force, arranged as in Fig. [TJ are of much smaller 
intensity than the MOT beams. Therefore, they are negligible when the MOT beams are ON. For the experimental 
detection of the synthetic Lorentz force, we turn the MOT beams OFF, as described in detail below. 

We need to measure the transverse force F y in dependence of the velocity v x of the atomic cloud. Thus, we must 
prepare a cloud with a given centre of mass (CM) velocity. For this purpose we use an additional pair of current coils 
that produce a bias magnetic field along the symmetry axis of the anti-Helmholtz coils, x in our notation here. The 
bias field moves the centre of the trap (the point where B = 0), which displaces the cloud approximately 1 mm along 
the x-axis. 

The measurement protocol is as follows, (i) We load the trap with the bias magnetic field on. (ii) At t = —2 ms we 
reverse the bias magnetic field by reversing the current, which suddenly shifts the centre of the trap. This introduces 
a force on the cloud due to the MOT beams. During the next 2 ms the cloud accelerates along the x-axis towards the 
new trap centre, (iii) Next we turn off the MOT cooling laser and all real magnetic fields. This moment corresponds 
to t = 0 in our presentation. The system is now simplified because Zeeman splitting of hyperfine levels is absent 
and the radiation force left, arising solely from the lasers implementing the synthetic Lorentz force, is not spatially 
dependent (it is only velocity-dependent), (iv) The cloud starts to expand because the trapping is absent, but it also 
moves in the xy plane due to both the initial velocity, v x {t = 0), and the radiation pressure force, (v) After some 
delay time t, the cooling laser is suddenly turned on and the cloud is imaged with the camera. From the trajectory 
traversed by the CM of the cloud (x(t),y(t)), we can find the force acting on the atoms. If we wish to image a cloud 
initially at rest we skip steps (i) and (ii). For a given delay time t, we repeat the measurement protocol 20 times in 
identical conditions, and subsequently average to obtain x(t) and y(t). The gravity is along the z-axis in our system; 
free fall of atoms due to gravity does not affect the motion in the xy plane. It should be stated that we perform these 
measurements first with 11) —>• 12) lasers OFF, and then with these lasers ON, keeping all other parameters identical. 
The difference in the path y(t), with lasers |1) |2) ON and OFF, gives us the transverse motion due to the synthetic 

Lorentz force. 

The experimental proof of the synthetic Lorentz force by radiation pressure. The results of the 
experiment are illustrated in Fig. [2] We show the trajectory of the cloud (x(t),y(t)), in the presence of the synthetic 
Lorentz force, for three initial velocities: v x = —0.3 m/s (squares), v x = 0 m/s (diamonds), and v x = 0.6 m/s (circles); 
v y = 0 at t = 0 in each run of the experiment. There is a difference in the magnitude of the initial v x for the positive 
and negative velocity [circles and squares in Fig. [51(a)], which is a result of our MOT retro-reflected geometry, and 
the way we accelerate the cloud in step (ii) of the protocol. In order to prepare a cloud with positive v x , the cloud 
is accelerated with the incoming MOT beam (coming from the laser side of the setup), whereas acceleration in the 
opposite direction is performed with the reflected beam which has smaller intensity. The reflected beam intensity is 
smaller due to the losses, which are a result of the passage of the incoming beam through the dense cloud (absorption), 
and partially due to reflection. Consequently, the negative initial velocity is smaller than the positive velocity. 

For the lasers implementing the synthetic Lorentz force, we use the following detuning values: for the first step at 
780 nm = — (5_>. = 2tt x 6 MHz, and for the second step at 776 nm 5-f = — 27 t x 3.5 MHz and = 27 t x 7.1 MHz. 
The two operational two-step pathways are P(-S—, f) and P(— *■,).), whereas the other two are far from resonance. 
The intensities of the beams used are Irso = 0.060 mW/cm 2 and I 776 = 2.9 mW/cm 2 , giving Rabi frequencies 
H 780 = 27r x 1.2 MHz and H 776 = 2n x 0.94 MHz. 

By inspection of Fig. [2] we see that the cloud travels along the x-axis by inertia, whereas it accelerates along 
the y-axis due to the synthetic Lorentz force. The direction of the force depends on the sign of the velocity v x 
(negative/positive v x gives positive/negative F y ), and the force is zero for a cloud at rest. We observe an asymmetry 
in the force F y , for the positive and negative velocity. In order to justify the exact choice of the detuning values and 
(5p and further investigate the observed asimmetry in the measured synthetic Lorentz force, we perform measurements 
in a slightly simplified configuration. 

Measurements in an auxiliary configuration. We inspect the force along the y direction arising from the two- 
step two-photon resonances, by using a configuration with three laser beams illustrated in Fig. [3] We block the beam 
pointing towards the negative y-direction, and measure F y arising from the remaining beam (the positive y-direction), 
which drives the |1) —> |2) transition. The force is measured for an atomic cloud with velocity zero, as a function of 
the detuning J-j-, see Fig. [3] Measurements are done for three different detuning values of the first-step 780 nm beams 
5 = — 5^. = 27t x 4,6, 8 MHz. The intensities of the lasers driving the transitions are now pso = 0.043 mW/cm 2 
and I 776 = 4.8 mW/cm 2 , giving Rabi frequencies H 780 = 27r x 1.1 MHz and H 776 = 27r x 1.2 MHz. The two maxima 
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FIG. 2. The trajectories of the CM of the atomic cloud in the presence of the synthetic Lorentz force, (a) x(t), and (b) y(t) 
for three different initial velocities, v x = 0.6 m/s > 0 (circles), v x = —0.3 m/s < 0 (squares), and v x = 0 m/s (diamonds); 
initial component of v y = 0 in all measurements. Accelerating motion along y is the signature of the transverse force F y , which 
depends on v x . The lines are fitted to the experimental data, linear fits for v x (t), and quadratic for v y (t). 


in Fig. [3] are profiles of the two-step two-photon resonances: the peak in the vicinity of 5-)- = — <5_>. > 0 corresponds to 
the excitation pathway P{— »,t), and the peak close to = —<5«_ < 0 corresponds to P(<— ,t)- 

Solid lines in Fig. [3] show the theoretically calculated profiles of F y (see Methods for details of the theoretical 
calculation). The agreement between theory and experiment is evident. All parameter values used in the theoretical 
calculation are taken from the experiment, except that the Rabi frequencies are reduced by 20%. This is reasonable 
because in the experiment, the absorption of laser beams across a dense atomic cloud reduces their intensity [2l| . 

It should be pointed out that the peaks in F y (S^) are slightly displaced from the values Sf = ±|<5_>.|, towards d-j- = 0. 
For example, for |5^| at 6 MHz, the maxima are at ±5.3 MHz. Moreover, the FWHM of the peaks is larger than the 
linewidth of the state |2). These two observations are a consequence of the laser linewidth, which is larger than the 
linewidth of the state |2) = |5H 5 / 2 ; F" = 4) (see Methods). The finite linewidth of |0) -A |1) lasers at 780 nm distorts 
the peaks as follows: the side of the peak closer to (5-j- = 0 is lifted up in comparison the opposite side of the peak. This 
follows from the fact that the two-step two-photon resonance is stronger when |(5_>.| is closer to zero, which is evident 
from Fig. [3] The exact positions of the peaks at ±5.3 MHz, explain the chosen detuning values for the second step at 
776 nm, which were used to obtain Fig. [2] <5q- = —27 t x 3.5 MHz and = 27r x 7.1 MHz. They are chosen such that 
the operational pathways are effectively equally detuned from the two-step two-photon resonance: P(— >-,4-) is detuned 
by 2 tt x (-5.3 ±7.1) MHZ= 2 t r x 1.8 MHZ, and P(<~A) is detuned by 2 tt x (±5.3 - 3.5) MHZ= 2t r x 1.8 MHZ. 


DISCUSSION 


Suppose that we repeat measurements corresponding to Fig. [3l but for an atomic cloud with mean velocity v x 
different from zero. The results of such measurements would be identical as for v x = 0, but the positions of the peaks 
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FIG. 3. Frequency scan of the two-step two-photon resonance in the auxiliary configuration. The calculated (solid lines) and 
measured force F y as a function of <5-f for the excitation with just three beams as shown. Measurements are performed for the 
cloud with initial velocity zero. The plots show resonances for three values of the detuning 5<_ = — = 2-7T x 4 MHz (circles), 

6 MHz (squares) and 8 MHz (diamonds). 


would correspond to the Doppler shifted detuning values S->. — kv x and 1 + kv x . Thus, because detuning can be 
mapped to velocity space, Fig. [3] can be reinterpreted as measurements for a fixed value of (5<_ = — and for three 
different velocities v x < 0 (4 MHz), v x = 0 (6 MHz), and v x > 0 (8 MHz). This is sketched in Figure 01 where we see 
that the two peaks separate (approach) each other for v x > 0 (v x < 0, respectively). 

We use Fig. [4] for a detailed explanation of the synthetic Lorentz force measured in Fig. [2] In measurements shown 
in Fig. H we have <5f < 0, which means that the positive force F y > 0 in Fig. [2] results from the left resonance peak 
in Fig. HJ Likewise, because <5^ > 0 was used for Fig. 0J the negative force F y < 0 results from the right resonance 
peak in Fig. HI The transverse force F y measured in Fig. [21 can be approximately identified with [FyiSf ) — F y (5i)\ Vx , 
as illustrated in Fig. 0] 

The choice of 5-j- and (h in Fig. [2] is such that the strength of the forces arising from the two peaks balance each 
other, giving F y = 0 for v x = 0 [Fig. 0Kb)]. Moreover, Sf and are on the slopes of the two peaks in Fig. 0Kb) for an 
atom with v x = 0, where the maximal magnitude of A F y /Av x is expected (detuning translates into velocity space via 
Doppler effect). For an atom with v x < 0 [see Fig. 0Ka)[, the two peaks approach each other, yielding greater force 
from the left peak, which results in F y > 0 for v x < 0, and the opposite for v x > 0 which yields F y < 0 as shown in 
Fig.0Kc). 

It should be pointed out that our scheme is inherently asymmetric. The intensity of the resonance peaks shown 
in Fig. 0] decreases when they separate (for v x > 0), in contrast to when they approach each other (for v x < 0). 
Thus, the net force along y is larger for negative v :r . than for the velocity of the same magnitude but with a positive 
sign. By using a scheme with four atomic levels 17|, one could remedy the asymmetry in the force, present in the 
three-level scheme. In addition, the initial density of the atomic cloud is not identical prior to acceleration in the 
positive/negative direction due to the slight asymmetry of the displaced trapping potential. This also affects the 
transverse force due to absorption and multiple scattering [21]. The inherent asymmetry in the force F y (v x ) and in 
the initial density of the cloud, is reflected in the asymmetry of the observed motion in Fig. [2] for v x < 0 and v x > 0. 
Importantly, all of these details only quantitatively affect the results, but not qualitatively. 


CONCLUSION 

The presented experiment, which demonstrates the synthetic Lorentz force by using radiation pressure, is performed 
in a classical cold atomic gas, prepared in a MOT, where the heating due to spontaneous emission does not prevent the 
observations. The concept of a force is natural in our classical laser cooled system [l8[. We would like to emphasize 
that our method is entirely different from the Berry phase approach, where the connection with the Lorentz force can 
be made in a semiclassical approximation [22| . 

In the outlook, we foresee many intriguing novel experiments based on the presented method. First, one could 
develop experiments using more sophisticated schemes, involving more atomic levels (see Ref. [l7|), to create a 
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FIG. 4. Interpretation of the synthetic Lorentz force from Figs. [Hand[2] via two-step two-photon resonances presented in Fig. 
[3] Sketch of the resonant peaks that would be obtained in the setup shown in Fig. [3] for three atomic velocities: (a) v x < 0, 
(b) v x = 0, and (c) v x > 0. Vertical dotted lines illustrate the values of the detuning used for the lasers aligned on the y -axis 
in Fig. |T| Arrows denote the direction of the force exerted by those lasers, and illustrate the way F y observed in Fig. [Harises 
as a function of velocity v x . See text for details. 


uniform synthetic magnetic field. Next, we plan to adjust our system for the observation of the predicted rotation of 
the cloud during expansion [l7) ■ One of the goals of this research is to build up a toroidal trap for cold atoms with a 
toroidal synthetic magnetic field, which holds potential to emulate the plasma in a tokamak. The proposed concept 
could be used for velocity selection of atomic beams, or for developing a novel kind of mass spectrometer for neutral 
atoms. We believe that our concept or an analogous scheme could be applicable in other systems, such as suspended 
nanoparticles. 


METHODS 

Experimental setup. The s 'Rb MOT is set up in the standard a + cr~ retro-reflected configuration, with beam 
diameters of 2 cm. The trap is vapour loaded in a glass cell which facilitates fast switching of the magnetic field. 
Cooling and repumper lasers are external cavity diode lasers (ECDL) delivering total powers of 80 mW and 20 mW, 
respectively. The cooling laser is typically 27 tx24 MHz (4T) red-detuned from the |5 Si/2 ;P = 2} —> \5P 3 /2;F' = 3) 
hyperfine transition. The repumper laser is in resonance with the | 5Si/ 2 ;F = 1} -4► \5Pi/2iF' = 2} hyperfine 
transition, thus keeping most of the population in the \5Si/ 2 ;F = 2) ground level. The quadrupole magnetic held 
gradient is 13 G/cm. The number of atoms in the trap is deduced by measuring the cloud fluorescence with a 
calibrated photodiode. 

For the implementation of the synthetic Lorentz force two additional ECDL lasers were introduced in the experiment: 
one at 780 nm driving |5 S 3 / 2 \F = 2) —> |5P 3 / 2 ;F' = 3) transition, and other at 776 nm driving |5P 3 / 2 ;P' = 3) —>■ 
|5P> 5 / 2 ; F" = 4) transition. Each laser beam is split into two beams which are sent on the rubidium cloud in counter- 
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propagating configurations as shown in Fig. [lj Frequency and intensity control is done separately for each of the four 
laser beams with acoustic-optical modulators (AOM). After being frequency shifted, we couple the beams to single 
mode polarization maintaining fibres, which ensures linear polarization (in the z-direction) and uniform intensity. 

All lasers used in the experiment are frequency stabilized by using modulation transfer spectroscopy. We modulate 
the laser diode current to modulate the laser frequency, which effectively increases the laser linewidth. This additional 
broadening is taken into account in theoretical calculations. The linewidths were checked by heterodyne beating of two 
stabilized lasers with similar locking parameters. For locking the laser tuned at 776 nm, we counter-propagate picked off 
beams from the 776 nm laser and the 780 nm laser through a heated 8 ' Rb glass cell, where the 780 nm laser populates 
the \5P 3 / 2 ;F' = 3) level. We observe an absorption signal resulting from the | 5 P 3 / 2 ;-F' = 3) —> \5D 5 / 2 -,F" = 4} 
transition. This signal is mixed with the modulation signal from the 780 nm laser to obtain a frequency locking error 
signal. Therefore, there is no need to modulate the 776 nm laser frequency to stabilize it. 

Theoretical calculation of the radiation pressure force. To calculate the force plotted in Fig. [3l we first 
solve the optical Bloch equations (OBEs) to find the density matrix p [Tij ]: 



dt 



nm 


r inmPnm ? 


n/m, 


dpnn 

dt 




E m >E n 


E m <E , 


r mnPnrn 


(i) 


where n,m= 0,1,2. The Hamiltonian H = 7F°) + TP 1 ) describes a three-level system interacting with laser fields. 
ff(°) is represented by a diagonal matrix with elements = 0 , = Hlooi, and H^ = fi(woi + wi 2 ), whereas the 

interaction Hamiltonian is, in the dipole and the rotating wave approximation 0 , 
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O i(—k x x—ujoit—5i-t) 
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and zero otherwise. The linewidths, Ti 2 = 27r x 6.1 MHz, r 23 = 27r x 0.66 MHz, and the transition frequencies, 
woi = k x c and u>i 2 = k y c , where k x = 27r/780 nm -1 and k v = 27t/ 776 nm -1 , correspond to the experiment; Ti 3 = 0; 
the coherences are 7 ij = \ (Fj, + Tjj ), where F,;, = Because transition |0) —> |1) is driven with two lasers 

of different detuning values <5<_ = — there is no stationary solution, and the level populations oscillate with 
the frequency |<5^| ~ 6 MHz (the phase of these oscillations depends on x, but this is not reflected in the force). 
The force on an atom is calculated by employing the Ehrenfest theorem [18|, and averaging over the oscillations: 
Fy{8^,8^,5^) = (—Tr\pdH/dy))f The average over time (■■■)$ is justified because the oscillations at frequency 
|<5_>| ~ 6 MHz are much faster than the characteristic scale for measuring the force, Tp, which is on the order of a few 
milliseconds. 

The outlined procedure for calculating the force assumes that the lasers are perfectly monochromatic. However, 
there is finite laser linewidth that should be taken into account to quantitatively describe the two-step two-photon 
resonances in Fig. [3] i.e., F y (S- (-)■ The spectral profile of the diode laser can be described by a Gaussian in the 
frequency domain, oc exp((w — wo) 2 /2cr 2 ). In the experiment, the modulation transfer spectroscopy for laser locking is 
used, and the central laser frequency is modulated with the amplitude d ~ 27r x 1 MHz and frequency 77 = 27r x 14 kHz: 
uio(t) = wq + dsinr/t. Thus, the laser spectral profile, over time-scales larger than I/ 77 , is given by 


p(w;wo) 


Poyd r o+d exp (—(a; - di 0 ) 2 /2cr 2 ) ^ 

<J\p2jF > Jw 0 —d |^f|t(wo) 

Po /'“ 0+d exp(-(w - w 0 ) 2 /2cr 2 ) 

- / / - , = CLUJ 0 

dv27r 3 Ju 0 -d A / ti 0 -Ld 0 \ 2 


(3) 


where po is the normalization factor. 

The detuning values of 780 nm lasers oscillate as follows: <5_>. = <5_> + d sin r/t and (5<_ = + dsinrjt. Note that 

r/ | A_j. |, and at the same time I /77 is sufficiently smaller than 77 -■ This means that we can employ the separation of 

scales to take into account finite laser linewidth. We calculate the resulting force F y {8 -j-) on the cloud by averaging 
F y {8 (5<_, <5-|-) with the appropriate laser profile distribution in the frequency domain: 

/ OO 

p(w;w 0 i + 8^)F' v (lo - ujqi,u> - w 0 i + 2|5_>|, 8-f)dw. 

-OO 


(4) 
















The width cr was determined to be 2n x 1.5 MHz by fitting to the experimental profiles. The Rabi frequencies used 
in the calculation are discussed in the Results Section. 
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